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Abstract — Temperature and doping dependencies of electron mobility in InP, InAs,GaP and GaAs structures have been 
calculated using an iterative technique. The following scattering mechanisms, i.e, impurity, polar optical phonon, 
acoustic phonon and piezoelectric are included in the calculation. The electron mobility decreases monotonically as the 
temperature increase from 100 K to 500 Kfor each material which is depended to their band structures characteristics. 
The low temperature value of electron mobility increases significantly with increasing doping concentration. The 
iterative results are in fair agreement with other recent calculations obtained using the relaxation-time approximation 
and experimental methods. 
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I. INTRODUCTION 

The problems of electron transport properties in semiconductors have been extensively investigated both 
theoretically and experimentally for many years. Many numerical methods available in the literature (Monte Carlo method, 
Iterative method, variation method, Relaxation time approximation, or Mattiessen's rule) have lead to approximate solutions 
to the Boltzmann transport equation 1 1-2]. In this paper iterative method is used to calculate the electron mobility of InAs, 
GaAs, InP and Gap. Because of high mobility InP has become an attractive material for electronic devices of superior 
performance among the Illphosphates semiconductors. Therefore study of the electron transport in group Illphosphates is 
necessary. GaP possesses an indirect band gap of 0.6 eV at room temperature whereas InP have a direct band gap about 1.8 
eV, respectively. InP and InAs offer the prospect of mobility comparable to GaAs and are increasingly being developed for 
the construction of optical switches and optoelectronic devices. 

The low-field electron mobility is one of the most important parameters that determine the performance of a field- 
effect transistor. The purpose of the present paper is to calculate electron mobility for various temperatures and ionized- 
impurity concentrations. The formulation itself applies only to the central U valley conduction band. We have also consider 
band nonparabolicity and the screening effects of free carriers on the scattering probabilities. All the relevant scattering 
mechanisms, including polar optic phonons, deformation potential, piezoelectric, acoustic phonons, and ionized impurity 
scattering. The Boltzmann equation is solved iteratively for our purpose, joint I in orpot tin th effi I oi ill the scattering 
mechanisms [3-4].This paper is organized as follows. Details of the iteration model , the electron scattering mechanism 
which have been used and the electron mobility calculations are presented in section II and the results of iterative 
calculations carried out on Inp,InAs, GaP, GaAs structures are interpreted in section III. 

II. MODEL DETAIL 

To calculate mobility, we have to solve the Boltzmann equation to get the modified probability distribution 
function under the action of a steady electric field. Here we have adopted the iterative technique for solving the Boltzmann 
transport equation. Under the action of a steady field, the Boltzmann equation for the distribution function can be written as 

^ + v .V r/+ ^.V, / = (^) fo „ a) 

dt h dt 

Where (df I ' dt) coll represents the change of distribution function due to the electron scattering. In the steady-state and 
under application of a uniform electric field the Boltzmann equation can be written as 

f.v,/ = (f)„„ 

h dt 

Consider electrons in an isotropic, non-parabolic conduction band whose equilibrium Fermi distribution function is fo(k) in 
the absence of electric field. Note the equilibrium distribution fo(k) is isotropic in k space but is perturbed when an electric 
field is applied. If the electric field is small, we can treat the change from the equilibrium distribution function as a 
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perturbation which is first order in the electric field. The distribution in the presence of a sufficiently small field can be 

written quite generally as 

/(*) = /„(*) + /,(*) cos 6 (3) 

Where is the angle between k and F and fj(k) is an isotropic function of k, which is proportional to the magnitude of the 

electric field. f(k) satisfies the Boltzmann equation 2 and it follows that: 

h dt 

In general there will be both elastic and inelastic scattering processes. For example impurit) scattering is elastic and 

acoustic and piezoelectric scattering arc clastic to a good approximation at room temperature. However, polar and non- 
polar optical phonon scattering arc inelastic. Labeling the elastic and inelastic scattering rates with subscripts el and inel 

respectively and recognizing that, for any process /'. s ell (A\ k) - s dl (/\, A") equation 4 can be written as 

^5 + 1 IV. W ISIJI - /„) + S^Qd'k' 



^ = i fcos P /;b:a-/ n )+s/ n rt-/ l j \s.(i-n+s.r <rt' 



. (5) 



: for elastic scattering. Equation 5 may be 



+ Xj/i'cos^[n-l][5L,(l-/o) + S,„,/ ]rf 3 r 



Note the first term in the denominator is simply the momentum relaxation r; 
solved iteratively by the relation 

-eF df 
f (#) = _JL 8k 

Xjd-cos* )s el d 3 k'+^[s m ji-f ) + s;, u ,X]d'k' 

Where / ln (k) is the perturbation to the distribution function after the n-fh iteration. It is interesting to note that if the initial 
distribution is chosen to be the equilibrium distribution, for which / j (k) is equal to zero, we get the relaxation time 
approximation result after the first iteration. We have found that convergence can normally be achieved after only a few 
iterations for small electric fields. Once / [ (k) has been evaluated to the required accuracy, it is possible to calculate 
quantities such as the drift mobility \x , which is given in terms of spherical coordinates by 



j(k 3 /l + 2aF)f l d 3 k 

]k 2 f d 3 k 



Here, we have calculated low field drift mobility in III-V structures using the ite: 
is electron-phonon and electron-impurity scattering mechanisms will be discussed. 



e technique. In the following 



Deformation potential scattering 

The acoustic modes modulate the inter atomic spacing. Consequently, the position of the conduction and valence 
band edges and the energy band gap will vary with position because of the sensitivity of the band structure to the lattice 
spacing. The energy change of a band edge due to this mechanism is defined by a deformation potential and the resultant 
scattering of carriers is called deformation potential scattering. The energy range involved in the case of scattering by 
acoustic phonons is from zero to 2vk Pi , where v is the velocity of sound, since momentum conservation restricts the 
change of phonon wave vector to between zero and 2k. where k is the electron wave vector. Typically, the average value of k 
is of the order of 10 7 cm" 1 and the velocity of sound in the medium is of the order of 10 5 cms" 1 . Hence, 2vkfo ~ 1 meV, 
which is small compared to the thermal energy at room temperature. Therefore, the deformation potential scattering by 
acoustic modes can be considered as an elastic process except at very low temperature. The deformation potential scattering 
rate with either phonon emission or absorption for an electron of energy £ in a non-parabolic band is given by Fermi's 
golden rule as [3,5] 



_ 42D 2 ac (m* 2 m*) ll2 K B T JE(l + aE) 



fa£) 2 +l/3(a£) 2 ] 



Where D ac is the acoustic deformation potential, p is the material density and a is the non-parabolicity coefficient. The 
formula clearly shows thai the acoustic scattering increases with temperature 

Piezoelectric scattering 

The second type of electron scattering by acoustic modes occurs when the displacements of the atoms create an 
electric field through the piezoelectric effect. The piezoelectric scattering rate for an electron of energy E in an isotropic, 
parabolic band has been discussed by Ridley [4] . 
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The expression for the scattering rate of an electron ii 
be written as [3,5]: 



-y 1/2 (l + 2a£J 



n-parabolic band structure retaining only the important terms c. 



of the material and K m , is the dimensionless so called average 



PZ l4inh 1 E s 

Where E s is the relative dielectric 
electromechanical coupling constant. 

Polar optical phonon scattering 

The dipolar electric field arising from the opposite displacement of the negatively and positively charged atoms 
provides a coupling between the electrons and the lattice which results in electron scattering. This type of scattering is called 
polar optical phonon scattering and at room temperature is generally the most important scattering mechanism for electrons 
in III-V .The scattering rate due to this process for an electron of energy E in an isotropic, non-parabolic band is [3-5] 



*„(')=- 



•I 2m a 



1 1 



\ + 2aE) 

'1M 



(10) 



8wr„/7 ^ 

f po (e,e'){n op ,n oi ,+i} 

Where E = E'±_w po is the final state energy phonon absorption (upper case) and emission (lower case) and Nop is 
the phonon occupation number and the upper and lower cases refer to absorption and emission, respectively. For small 
electric fields, the phonon population will be very close to equilibrium so that the average number of phonons is given by the 
Bose- Einstein distribution. 

Impurity scattering 

This scattering process arises as a result of the presence of impurities in a semiconductor. The substitution of an 
impurity atom on a lattice site will perturb the periodic crystal potential and result in scattering of an electron. Since the mass 
of the impurity greatly exceeds that of an electron and the impurity is bonded to neighboring atoms, this scattering is very 
close to being elastic. Ionized impurity scattering is dominant at low temperatures because, as the thermal velocity of the 
electrons decreases, the effect of long-range Coulombic interactions on their motion is increased. The electron scattering by 
ionized impurity centers has been discussed by Brooks Herring [6] who included the modification of the Coulomb potential 
due to free carrier screening. The screened Coulomb potential is written as 



V(r) = 



exp(-q r) 
4ae e s r 
Where £ s is the relative dielectric 
degenerate conditions is given by 



of the material and q is the ii 



g length, which under n 



i; KJ 



Where n is the electron density. The scattering n 
N/(l + 2aE) \ T „ n ^_bA 



i for an isotropic, nonparabolic band 



h'q 
Where N t is the impurity 

III. RESULTS 

We have performed a series of low-field electron mobility calculations for GaP, InP, GaAs and InAs materials. 
Low-field motilities have been derived using iteration mefhode.The electron mobility is a function of temperature and 
electron concentration . 



Figures 1 show the comparison the electron mobility depends on the Temperature at the different electron 
n bulk Gap.GaAs,InP and InAs materials, figures 1 show that electrons mobility at the definite temperature 
300 k for the InAs semiconductors is gained about 8117cm2v-ls-l and for GaAs, InP, Gap about 4883,4702,760 cm2v-ls- 
l.also electrons mobility decrease quickly by temperature increasing from 200k to 500 k for all the different electron 
s because temperature increasing causes increase of phonons energy too. So it causes a strong interaction 

58" 
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between electrons and these phonons that its result is increase of electrons scattering rale and finally decrease of electrons 
mobility. Our calculation results show that the electron mobility InAs is more than 3other materials this increasing is because 
of small effective mass. 




Fig 1. Changes the electron mobility Function 
electron 



of temperature in bulk InAs, GaAs, InP and GaP at the 
1016 (cm-3). 

Figure 2 shows comparison the electron mobility of InAs, GaAs, InP and GaP at the different Temperature .Our 
calculation results show that the electron mobility InAs is more than GaAsJnP and GaP 
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Fig 2. Changes the electron mobility Function in terms of room temperature in bulk InAs, GaAs, InP and GaP at the 
different electron concentration. 

Figures 3 shows comparison the electron mobility depends on the electron concentration at the different 
Temperature in bulk GaP and InP materials. Semiconductors mobility decrease by electrons concentrations increasing 
because elections increasing causes increase of ionized impurity centers in crystals that it causes times more electrons under 
the influence of the Coulomb potential of impurity centers located that its result is increase of electrons scattering rate and 
finally decrease of electrons mobility 
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Fig 3. Changes the electron mobility function in terms of electron 
Temperature 
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Fig 4. Changes the electron mobility Function it 



is of temperature in bulk GaAs and InAs at the different electron 



Figures 4 shows comparison the electron mobility depends on temperature at the different electron concentration in bulk 
InAs and GaAs materials. .Our calculation results show that the electron mobility InAs is more than GaAS 

IV. CONCLUSION 

1. InAs semiconductor having high mobility of GaAs,GaP and InP because the effective mass is small compared 
with all of them 

2. The ionized impurity scattering in all the semiconductor InAs,GaAs,InP and GaP at all temperatures is an 
important factor in reducing the mobility. 
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